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In extension of our previous study of KNLO3 by the semiempirical Hartree-Fock method we present 
parameterization and total-energy results for non-ferroelectric KTaC>3 as a pure crystal (concentrat- 
ing on the frozen phonon calculations) and that with Li impurities. The magnitudes of off-center 
Li displacements and the relaxation energies related to re-orientation of Li are calculated and com- 
pared with experimental estimates and earlier calculation results. The spatial extent of lattice 
relaxation around Li impurities and contributions from different neighbors to the relaxation energy 
are discussed. 
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I. INTRODUCTION 

Potassium tantalate is a so-called incipient ferroelectric 
which develops a high value of the dielectric constant and 
considerable softening of a zone-center transverse-optical 
(TO) phonon, typical to other ferroelectric materials, but 
stays in a paraelectric state even at the lowest tempera- 
tures. This paraelectric state can be modified by doping, 
that may (in case of Nb substituting Ta) bring the sys- 
tem to a true ferroelectric state, or (for doping with a 
sufficiently small amount of Li or Na that enter the K 
sublattice) gives rise to a more complicated state, that is 
often referred to as orientational glass (see, e.g., Refs. E. 
^ for a review) . As was discovered by Yacoby and JustB, 
a Li ion substituting K in KTa03, because of its smaller 
ionic radius, gets spontaneously displaced along one of 
the (six) [100] directions, towards an octahedral intersti- 
tial. This effect has been further studied in Refs. [|, ||, 
^. An experimental estimate of the displacement magni- 
tude possible from 7 Li NMR measurements (based on the 
measured values of the local electric field gradient, inter- 
preted in terms of ionic displacement) is 0.86 A according 
to Bcwsa et alu. Another estimate cited by van der Klink 
et al. - , in addition to the former one, is -1.2 A. Hochli 
et al.U report the displacement mag nitude of 1.1 ±0.1 A. 
These values are therefore to some extent dependent on 
the model in question that provides the fit to the data 
directly measurable by NMR.p 

Van der Klink and KhannaQ calculated the energetics 
of the [100] Li displacement in a polarizable point-charge 
model, with Coulomb and polarization energy treated as 
those from point charges (of nominal ionic values), with a 
selected set of polarizabilities. The equilibrium displace- 
ment (1.35 A) agrees well with the NMR-based estima- 
tion of —1.2 A, but the model is probably too crude for 
handling a delicate balance between long-range electro- 
static forces and short-range effects of chemical bonding 
in an incipient ferroelectric system. Particularly, it be- 
came well known since then that effective charges in per- 
ovskites deviate considerably from nominal ionic values, 
and the lattice relaxation for instance in the KTaOsiLi 
system gives a substantial contribution to the energy low- 
ering on a Li displacement. 



Several studies on Li in KTa03 have been done with 
the use of the shell model which is more sophisticated 
than the models dealing with^oint charges or point 
dipoles. Stachiotti and MigoniQ calculated equilibrium 
displacements of Li and neighboring atoms making use of 
the static lattice Green function method formulated on 
top of the nonlinear shell model, with anisotropic core- 
shell couplings of the oxygen ion. The drawback of this 
approach is that the anisotropic (fourth-order) interac- 
tion was allowed for the O-Ta bond only, but otherwise 
the model was fully harmonic. This seems to be hardly 
satisfactory for quantitative estimates, taking into ac- 
count the large magnitude of the Li displacement. More- 
over, there was an ambiguity in fitting the Li-0 potential 
of the Born-Mayer tyge, based either on a study of lattice 
dynamics in LiKSO^, or on the studw-pf defect energies 
of Li+ substituting Ba 2+ in BaTioJil The results of 
Ref. |[ obtained with these two sets of parameters, dif- 
fer considerably: the first set leads to a non-displaced 
position of Li, whereas the second set provides the [100] 
displacement bx 1.44 A that seems to be too large. 

Exner et aZ.Eil performed another simulation of doped 
KTaO,3 within the shell model, with the polarizability 
of the oxygen ion treated as isotropic, but, on the other 
hand, with anharmonic effects accounted for. Also the 
choice of short-range interaction was different in Ref. |ll| 
from that of Ref. ||. The equilibrium off-center Li dis- 
placement was found to be ~0.64 A, and the relaxation 
pattern of near neighbors to Li impurity essentially dif- 
ferent from that provided by Staciotti et-alB. 

First-principles supercell calculationsEj by the full- 
potential linear muffin-tin orbitals method (FP-LMTO) 
made it possible to study the total energy as function of 
off-center Li displacements from its K-substituting posi- 
tion in KTaOa. This calculation, although not depen- 
dent on the choice of the particular interaction model or 
the fitting of the potential, had nevertheless its technical 
limitations. For one thing, the computational effort rises 
quite fast with the supercell size. The largest supercell 
treated in Ref. [l^ included 2x2x2 primitive cells, i.e. 
only 40 atoms. Another difficulty relates to the muffin- 
tin geometry used in FP-LMTO calculations (see, e.g., 
Refs. 0, 04). This assumes all atoms to be circumscribed 
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by non-overlapping spheres of reasonable size, so that the 
potential and charge density are expanded in spherical 
harmonics inside such spheres, and the interstitial region 
is treated in some different way. The large magnitude of 
the Li off-center displacement and, as a result, of lattice 
relaxation makes the universal non-overlapping packing 
for all geometries possible only with quite small spheres, 
that is disadvantageous for the numerical accuracy. 

The method of our choice in the present study is 
therefore a tight-binding scheme with the basis of atom- 
centered orbitals, that makes no use of muffin-tin geom- 
etry whatsoever. At the same time, we want the scheme 
to be not much computationally demanding, in order to 
keep the calculations for large supercells feasible. What 
in some aspects closes the gap between empirical models 
(applicable to very large systems) and accurate ab ini- 
tio schemes (that become very computationally demand- 
ing unless applied to supercells of quite modest size) are 
semiempirical methods which are able to produce reli- 
able quantitative predictions, based on a limited num- 
ber of basic (system-dependent) tunable parameters. We 
use the Intermediate^feglect of the Differential Overlap 
(INDO) methodtallaO, which has been applied very suc- 
cessfully for the study of defects, both r 
on the surface, in many oj 
as well as semiconductor: 



materials 

e recently applied this 
method to the study of KNbOaE3. The method is essen- 
tially a simplified implementation of the Hartree-Fock 
formalism, with its obvious drawback of underestimat- 
ing the correlation effects. The latter seems however to 
be of little importance when one studies an ionic insula- 
tor as KTaC>3 and not, say, a metallic system. Moreover, 
the appropriate choice of parameters makes it possible to 
obtain reasonable trends in the total energy as function 
of displacements even if the band structure is somehow 
distorted in the Hartree-Fock calculation (see Ref. for 
the discussion to this point for KNbOs). In tuning the 
INDO parameters, we primarily use earlier FP-LMTO 
calculations as a benchmark, also with some reference to 
experimental data (e.g., for phonons in KTaOs). 

The objectives of the present study are the total energy 
as a function of off-center Li displacements in KTaOs, the 
90°-energy barrier in the hopping motion of displaced Li, 
and the lattice relaxation around the Li impurity. As a 
by-product, we provide the description of pure KTa03 
by the INDO method, including the calculation of L 
phonons. The paper is organized as follows. In Sec. |J, 
we discuss the choice of INDO parameters for our cal- 
culations of pure and Li-doped KTa03. In Sec. OlL we 




frequencies and eigenvectors in KTa03. In Sec. 
discuss the off-center Li displacement and the la' 
laxation around a Li impurity in KTa03. 
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II. INDO PARAMETERIZATION FOR KTaO a 

The description of the INDO approximation to the 
Hartree-Fock-Roothaan method is given in Refs. ||, [IT]. 
The major formulae defining the parameters of the calcu- 
lation are also cited in Ref. g5|, where the choice of param- 
eters relevant for the calculation of KNb03 is discussed. 
As in Ref. |25|, we performed the present calculations with 
the CLUSTERD-cpmputer code by Shidlovskaya, Shlugcr 
and StefanovicrfiaO. The parameters obtained there for 
K and O atoms are retained in the present calculation. 
In the following, we refer to Ref. ^5| for the description 
of relevant parameters. For KTaOs, we had to specify 
£ M (parameter of the Slater exponent), E neg (central en- 
ergy position), /3 M (parameter of the resonance interac- 
tion with other states) and (occupation number) for 
the Ta and Li states. Moreover, the two-center param- 
eter olab had to be specified for Ta and Li in combina- 
tions with_other atoms. r-The INDO calculations for Li 
halidesE3E3 and Li2Si03E3 have been done earlier, and 
some parameters tabulated. In contrast to these earlier 
works, we preferred to use both Li 2s and 2p states in 
the basis set, since the bonding in KTa03:Li is much 
less ionic than in Li halides. Moreover, a very delicate 
balance of the interactions resulting in a ferroelectric in- 
stability favors the use of a more extended basis set. The 
recommended values of one-center Li parameters—which 
are provided in the description of the INDO codeO were 
found to be a reasonable choice (see below), and some 
discussion as for the choice of two-center ao-Li parame- 
ter is given in Sec. Evl. As for the Ta-related parameters, 
we took E neg and P^ values (with the exception of E neg 
for Ta 5d) the same as for corresponding (5s, bp and 4c?) 
Nb states in KNb03, based on the similarity of the band 
structures. {Ab initio calculations of the band structure 
of both compounds have been performed earlier, e.g., by 
Neumann et alr3). 

The final choice of the INDO parameters for KTa03 
was done based on a calculation of frozen L-TO phonons. 
Since the three phonon modes which belong to the T\ u ir- 
reducible representation are different combinations of the 
vibrations of K, Ta, and O, they provide a good indica- 
tion of the accuracy with which all essential interatomic 
interactions are treated. The TO phonon frequencies are 
measured in a number of experiments, and their eigen- 
vectors determined in ab initio calculations. This makes 
the adjustment of aAB and /3 M parameters, which affect 

TABLE I. One-center INDO parameters for Ta and Li 
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Orbital 


C (a.u.- 1 ) 


£ ncg (eV) 
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Ta 6s 


2.05 


0.0 


20.0 


0.1 




Ta 6p 


2.05 


-2.0 


20.0 


0.0 




Ta 5d 


1.70 


23.50 


16.0 


0.6 




Li 2s 


1.10 


5.0 


0.6 


0.1 




Li 2p 


0.80 


0.9 


0.6 


0.02 
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TABLE II. Calculated T-TO frequencies and eigenvectors of the T\ u modes in cubic KTaC>3. 
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u> calc. (cm ) 


uj exp. (cm ) 


K 


Ta 




^2,3 






0.06 


-0.46 


0.71 


0.37 


86 a 


25-106 c ; 81 d ; 85 e 


0.11 


-0.50 


0.47 


0.51 


80 b 




-0.92 


0.31 


0.14 


0.13 


202 a 


196-199 c ; 199 d ; 198 e 


-0.91 


0.28 


0.14 


0.17 


172 6 




-0.01 


0.13 


-0.65 


0.53 


500 a 


551-550 c ; 546 d ; 556 e 


0.01 


0.09 


-0.83 


0.39 


528 6 





a Present work. 

b Full-potential LAPW calculation of Ref. 

c Infrared reflectivity measurements at 12—463 K, Ref. B6[ 

d Hyper-Raman scattering measurements at room temperature, Ref. [33]. 

e Raman scattering measurements at room temperature (soft mode) and at 10 K, Ref. 



the shape of the total energy hypersurface as function 
of a particular displacement and changes the frequencies 
and eigenvectors, a relatively easy task. To begin with, 
we took as a benchmark the total energy difference as 
function of individual displacements, determined in the 
full-potential LMTO calculations However, finally we 
aimed at obtaining the phonon eigenvectors with appar- 
ently even better accuracy, as calculated in Ref. |3(] by 
the full-potential linear augmented plane waves (LAPW) 
method. The values of the INDO parameters we found 
to be best suited for the study of KTaOsiLi (in addition 
to those published already in Ref. |2^) are listed in Ta- 
ble |. The two-center parameters olab which account for 
the non-point character of the interaction of a valence or- 
bital at the atom A with the core of atom B (see explicit 
definition in Refs. fj]], ||) are 0.023, 0.15, 0.39 and 0.58 
for A—O and B=L\, O, K and Ta, correspondingly, and 
zero for ^4=Ta, K, or Li. 

The band gap in KTaC>3 as calculated by the INDO 
method is 6.7 eV, that is close to, but larger than, 6-J 
eV as obtained in a similar calculation for KNbOsEi 
The absolute value of the gap in the one-electron energy 
spectrum is known to come out systematically larger in 
the Hartree-Fock formalism as compared with spectro- 
scopic data, because the unscreened Coulomb interaction 
shifts the unoccupied states too high in energy. Never- 
theless, the difference between the gap values for two 
compounds is in agreement both with the experimental 
estimates (3.3 eV for KNb0 3 vs. 3.8 eV for KTa0 3 ) and 
the results of the calculations donaia the local density 
approximation (1.4 eV vs. 2.1 eVOEl). The (static) ef- 
fective charges for KTa03, as estimated from the INDO 
calculations based on the Mullikan population analyzis, 
are +0.62 (K), +2.23 (Ta) and -0.95(0). This reveals 
a somehow increased ionicity for KTa03, in comparison 
with KNb03t3, in geperal agreement with the tendency 
pointed out by SinghtH 



III. T-TO PHONONS IN KTaO s 

T-TO phonons in the cubic perovskite structure are 
split by symmetry into three triple degenerate T\ u modes 
and one triple degenerate T2„ mode. Ab initio calcula- 
tions of frequencies and eigenvectors have been done by 
Postnikov et a/.tll and by SinghEx Whereas calculated 
frequencies are in reasonable agreement between both 
calculations and with the experimental data available, 
there is some disagreement in the estimations of calcu- 
lated eigenvectors, especially for the soft mode, as was 
discussed at length in Ref. [30[ The full-potential LAPW 
method is able to provide ultimately better accuracy than 
LMTO, due to a more extended basis set. Moreover, the 
indirect experimental indications of the displacements 
within the soft mode of KTa03&3 seem to be in agreement 
with the LAPW results. Therefore we aimed at reproduc- 
ing the eigenvectors of Ref. ^ with the proper choice of 
our INDO parameters. We performed our frozen phonon 
calculations at the experimental lattice constant (extrap- 
olated to zero temperature) of 3.983A. The elements of 
the dynamical matrix were found from the polynomial fit 
of the total energy as function of various displacements 
within the T± u mode. Calculated frequencies and eigen- 
vectors are shown in Table |j Oi stands for the oxygen 
atom at (^^0), and 02,3 for two (equivalent) atoms at 
(0|i) and (^0^), for the vibrations along [001]. K is at 
(000) and Ta at {\\\) of the cubic perovskite cell. ._. 

The agreement with the eigenvectors of SinghEJ is 
good, especially in what regards the relative displace- 
ment of potassium and tantalum with respect to each 
other and to the averaged displacement of the oxygen 
sublattice. The difference in the displacement patterns 
of Oi and 02,3 atoms within the soft mode seems, how- 
ever, to be slightly overestimated in our case. 

For the Ti u mode, our calculation gives the freauency 
of 260 cm -1 , as compared to 264 cm" 1 by SinghEJ and 
the experimental estimations of 264 cm" 1 (Ref. |32|) to 
274 cm -1 (Ref. |33|). This mode involves only the stretch- 
ing within the O2.3 sublattice, therefore the eigenvector 
is uniquely defined. The excellent agreement of the fre- 
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FIG. 1. [100] and [110] Li off-center displacements and the 
relaxation pattern of neighboring oxygen atoms. 



quency with the experimental value(s) suggests that our 
choice of the INDO parameters, that describe the 0-0 
interaction as mediated by the central Ta atom, is suffi- 
ciently good. 



IV. Li OFF-CENTER DISPLACEMENT AND 
LATTICE RELAXATION 

With the parameters of the INDO method properly 
tuned, one can obtain better accuracy in describing en- 
ergy trends and lattice relaxation in Li-doped KTa03 

jle point- 
On the 



than was possible with the simple po, 
charge modelQ or within the shell modeltlt. 
other hand, the INDO calculation allows to overcome the. 
problems of previous first-principles FP-LMTO studiestS 
in what regards small supcrcell size and the problems of 
the muffin-tin-spheres packing. Since the muffin-tin ge- 
ometry is not used in a tight-binding INDO scheme, the 
problems of spheres packing do not exist there; at the 
same time, the method is much less computationally de- 
manding and allows to treat larger supercells (all results 
discussed below refer to the 3x3x3 KTa03 supercell, 
i.e. that with 135 atoms in total, with one substitutional 
Li impurity). It is also possible to search for the opti- 
mized lattice distortion around the displaced impurity in 
the course of INDO calculation, that was not done in 
Ref. ||. 

In addition to the parameters which relate to the 
KTaO,3 bulk and which were tuned on the basis of frozen 
phonon calculations of Sec. [n], and the one-site Li pa- 
rameters borrowed from eadier INDO calculations for 
Li-containing ionic crystalalj'O, we had to specify the 
two-center parameter aQ_Lj, that enters the empirical 
expression for the interaction of the O valence electrons 
with the Li core and that effectively accounts for the 
finite size of the Li core and for the diffuseness of the O- 
related valence band states (see definition of this param- 
eter and the discussion in Ref. |l7|) . Whereas already the 
first choice of «O-Li = as proposed in Ref. [ij] provides 
a qualitatively correct effect on the equilibrium geome- 
try (Li displaces off-center along [100] by ^0.4 A), we 
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[110] Li displacement (A) [100] 

FIG. 2. Total energy gain as function of [100] and [110] Li 
displacements calculated by FP-LMTO method (Ref. fL3) and 
by INDO for the 3x3x3 supercell. The curves are shifted so 
as to match for the 0.3 A displacement. 



found that a fine tuning of this parameter has a consid- 
erable effect on the magnitude of the displacement, and 
especially on the energy gain associated to it. For in- 
stance, «o-Li— 0-02 sets the energy gain at ~60 meV at 
the ~0.6 A Li displacement, whereas ao-Li=0-04 low- 
ers the total energy by ~100 meV at the ~0.7 A dis- 
placement -i^The best fit to the results of the FP-LMTO 
calculation^ (that were done only for the Li displace- 
ment >0.3 A) occurs at ao-Li=0-023. We attempted 
to reproduce the FP-LMTO results for [100] and [110] 
Li off-center displacements - the former being related to 
the true ground-state configuration, and the latter to the 
saddle point between two adjacent displaced Li positions. 
The displaced Li ion with its neighboring oxygen atoms is 
shown schematically in Fig. |l|. The total energy as func- 
tion of [100] and [110] Li off-center displacements from 
FP-LMTO and from our present calculations is shown in 
Fig. | 

Another sensitive question is the dependence of the re- 
sults on the supercell size. The INDO method in the 
implementation we use produces in each iteration the 
one-electron energy spectrum in the T point only. The 
calculation is normally being done with an extended su- 
percell, so that the band dispersion over the correspond- 
ingly shrinked Brillouin zone is neglected. The enlarge- 
ment of a supercell in the calculation for a perfect crystal 
is in a sense analogous to choosing a finer mesh for the 
k-space integration in a conventional band structure cal- 
culation. When treating an impurity system, another 
aspect becomes essential as well, that is, the effect of in- 
teraction between impurities situated in adjacent super- 
cells should be kept negligible, or at least controllable. 
For the KTaOsiLi system, it means that the polariza- 
tion cloud associated with a single off-center displaced 
Li impurity should be, in the ideal case, fully within the 
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supercell chosen. The size of the polarized region associ- 
ated with the [100]-displaced Li ion was estimated, in the 
shell model calculation by Stachiotti and MigoniB to be 
about 5 lattice constants along the direction of displace- 
ment, with ~99% of the 'effective dipole' polarization 
being confined to nearest Ta-0 chains, that go parallel 
to the displacement. As is discussed below, the mag- 
nitudes of the atomic displacements and polarization in 
our present calculation is considerably smaller than those 
found in Ref. ||, and the relaxed neighbors to the Li im- 
purity are well within the 3x3x3 supercell. In order to be 
on safer side, we performed as well the calculations for a 
supercell doubled in the direction of Li displacement, i.e., 
6x3x3, with a single [100]-displaced Li atom. The equi- 
librium displacement in this case is 0.62 A, exactly as for 
the 3x3x3 supercell (see Fig. |^), with the energy lowering 
57.2 meV. The difference from the result for a 3x3x3 
supercell (62.0 meV) roughly represents the uncertainty 
related to the supercell size in our calculations. 

The magnitude of the Li off-center displacement najtn- 
rally agrees well with the FP-LMTO data (0.61 AH), 
since the latter was an important benchmark in our 
choice of INDO parameters. We failed however to obtain 
the ideal matching with the FP-LMTO results in what 
regards both the off-center displacement and the energy 
gain in both [100] and [110] directions, as is seen in Fig. g. 
Our equilibrium [100] off-center displacement is smaller 
than the estimate (1.35 A) of the polarizable point-charge 
modeH, or 1.44 A as calculated by Stachiotti and Migoni 
within the shell modc-13. On the other hand, our value 
is in good agreement with a more recent, and apparently 
more elaborately parameterized, shell model calculation 
by Exner et al. (0.64 A, Ref. |ll|). 

The energy gain due to the Li off-center displacement 
is not directly measurable in an experiment, but there 
are estimations for the 90°-energy barrier between, say, 
[100] and [010] -displaced positions to be 86 meVtl The 
hoppings between such adjacent positions may only oc- 
cur via the [110] saddle point (see Fig. ||). Our estimate 
of the energy difference between [100] and [110] minima 
is ^30.2 meV, roughly two times larger than in the FP- 
LMTO calculationEJ, but much less than the experimen- 
tal estimate. The origin of this discrepancy, as has been 
mentioned in Ref. [l^, is most probably related to the lat- 
tice relaxation around the displaced Li ion, that makes 
the net energy gain from the displacement larger, and 
the 90°-activation energy (involving now the displace- 
ment of many atoms) correspondingly higher. Indeed, 
the second harmonic generation-based estimates of the 
activation barriered reveal two types of processes, appar- 
ently one involving the lattice relaxation (with the barrier 
height 86.2 meV) and another one that is too fast for the 
lattice to follow, with the barrier 14.7 meV. 

In order to clarify this point, we performed a lattice 
relaxation of several shells of neighbors to the displaced 
Li ion, for the cases of [100] and [110] displacements. 
The relaxed coordinates of atoms are given in Tables III 
and IV, where the oxygen atoms are numbered consis- 
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FIG. 3. Total energy gain as function of [100] and [110] 
Li displacements without lattice relaxation (dashed line with 
open circles), and the total energy values after including the 
relaxation of three groups of nearest oxygen atoms. 



tently with Fig. [I]. The total energy values resulting from 
the gradual inclusion of neighbor relaxation are shown in 
Fig. |^. We found the relaxation of twelve nearest oxy- 
gen atoms essential, and the effect of relaxing nearest Ta 
and more distant atoms to be negligible, in what regards 
the effect on the total energy. The energy gain in the 
fully relaxed [100]-displaccd configuration, with respect 
to a non-relaxed central Li position, is 158.9 meV; the en- 
ergy gain in the relaxed [110]-configuration is 102.3 meV. 
Therefore, the enhancement of the excitation barrier due 
to relaxation effects is by a factor of two, but still not suf- 
ficient to reach experimentally expected ~86 meV. This 
discrepancy may be due to the fact that in reality the 90°- 
reorientation process of the impurity does not necessar- 
ily occur via the fully relaxed saddle-point configuration. 
Depending on the actual degree of relaxation around the 
saddle-point Li position, the barrier height is expected 
from Fig. || to be between ~57 meV (full relaxation at 
the saddle point) to ^127 meV (no relaxation). 

It should be noted that the spatial range of polarization 
is not so large in our calculation as follows from the shell 
model results of both Ref. || and [H], and the displace- 
ments of twelve nearest oxygen atoms around the Li ion 
are, generally, smaller in our case. There are also some 
qualitative differences in the displacement pattern. The 
largest striking, the 0(2) atoms (as labeled in Table III) 
follow the Li displacement in our calculation, whereas 
they move away from the displaced Li ion in both shell 
model calculations cited. This is probably due to an over- 
simplified parameterization of the interaction potential 
used in the shell model. In the INDO method, no special 
approximation is introduced for the description of the 
chemical bonding, therefore our results seem to be more 
reliable. The numerical values of the displacements may, 
however, be somehow refined in subsequent calculations 
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TABLE III. Relaxed atomic positions for the [100] Li displacement as calculated by INDO. 



Atom Lattice coordinates Displacement 

Li A~ x A x = 0.1550 

4x0(1) l+Az | + A x = -0.0045; A y = -0.0105 

4xO(2) A x \ + &y \ + ^v Aa; = 0.0070; A y = —0.0026 

4xO(3) \ + ^v A * = -0.0020; A y = 0.0020 



TABLE IV. Relaxed atomic positions for the [110] Li displacement as calculated by INDO. 



Atom 



Lattice coordinates 



Displacement 



Li 

lxO(l) 
4xO(2) 
2xO(3) 
4xO(4) 
1x0(5) 



A x 
\ + A x 
\ + A x 

\+A x 

-\+A x 
-|+A X 



Ax 
\+A x 

Ay 

-\+^y 



-\+A x 








-i+A z 



Az = 0.0760 
A x = -0.0090 
A x = -0.0060; A y = 0.0030; A z = -0.0080 



-0.0020; A„ 



A x = 0.0003, A„ 



: 0.0001 

: ~0 



0.0060 

A ^ = 0.0003 



with larger unit cells. 



V. SUMMARY 

As an extension of our previous study of ferroelectric 
KNbC>3 with a semiempirical INDO method, we per- 
formed calculations for pure paraelectric KTaC>3, con- 
centrating on TO phonon frequencies as a benchmark 
for fine tuning of our INDO parameterization. In a sc- 
ries of supercell calculations for Li-doped KTa03, tuned 
in such a way as to reproduce the energetics of the Li 
off-center displacement previously found in FP-LMTO 
calculations, we analyze the relaxation of near neighbors 
to the Li impurity, and the impact of this relaxation on 
the reorientational energy barriers. The relaxation pat- 
tern in some aspects differs from that calculated earlier 
within the shell model. The study of the interaction be- 
tween Li impurities in a polarized lattice seems feasible 
with the method used. 
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